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Abstract 

A search for diphoton events with large missing transverse momentum has been performed us- 
ing proton-proton collision data at yfs = 7 TeV recorded with the ATLAS detector, corresponding to 
an integrated luminosity of 4.8 fb 1 . No excess of events was observed above the Standard Model 
prediction and model-dependent 95% confidence level exclusion limits are set. In the context of a 
generalised model of gauge-mediated supersymmetry breaking with a bino-like lightest neutralino of 
mass above 50 GeV, gluinos (squarks) below 1.07TeV (0.87 TeV) are excluded, while a breaking 
scale A below 196 TeV is excluded for a minimal model of gauge-mediated supersymmetry breaking. 
For a specific model with one universal extra dimension, compactification scales 1/R < 1 .40 TeV are 
excluded. These limits provide the most stringent tests of these models to date. 
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Abstract 

A search for diphoton events with large missing transverse momentum has been performed using proton-proton collision 
data at y/s = 7TeV recorded with the ATLAS detector, corresponding to an integrated luminosity of 4.8fb~ 1 . No excess 
of events was observed above the Standard Model prediction and model-dependent 95 % confidence level exclusion limits 
are set. In the context of a generalised model of gauge- mediated supersymmetry breaking with a bino-like lightest 
neutralino of mass above 50 GeV, gluinos (squarks) below 1.07 TeV (0.87 TeV) are excluded, while a breaking scale A 
below 196 TeV is excluded for a minimal model of gauge-mediated supersymmetry breaking. For a specific model with 
one universal extra dimension, compactification scales 1/R < 1.40 TeV are excluded. These limits provide the most 
stringent tests of these models to date. 



1. Introduction 

This Letter reports on a search for diphoton (77) 
events with large missing transverse momentum (i?™ lss ) 
in 4.8 fb -1 of proton-proton (pp) collision data at y/s = 
7 TeV recorded with the ATLAS detector at the Large 
Hadron Collider (LHC) in 2011, extending and supersed- 
ing a prior study performed with 1 fb _1 [l[. The results are 
interpreted in the context of three models of new physics: 
a general model of gauge-mediated supersymmetry break- 
ing (GGM) pHl, a minimal model of gauge-mediated su- 
persymmetry breaking (SPS8) and a model with one 
universal extra dimension (UED) (B-Q ■ 

2. Supersymmetry 

Supersymmetry (SUSY) 0-[l7| introduces a symmetry 
between fermions and bosons, resulting in a SUSY part- 
ner (sparticle) with identical quantum numbers except a 
difference by half a unit of spin for each Standard Model 
(SM) particle. As none of these sparticles have been ob- 
served, SUSY must be a broken symmetry if realised in na- 
ture. Assuming i?-parity conservation [18-22], sparticles 
are produced in pairs. These would then decay through 
cascades involving other sparticles until the lightest SUSY 
particle (LSP), which is stable, is produced. 

In gauge-mediated SUSY breaking (GMSB) models [22 
the LSP is the gravitino G. GMSB experimental sig- 
natures are largely determined by the nature of the next- 
to- lightest SUSY particle (NLSP). In this study the NLSP 
is assumed to be the lightest neutralino Xi- F° r studies 
with the lightest stau as NLSP, the reader is referred to 
Refs. [29|, [3(|. Should the lightest neutralino be a bino 
(the SUSY partner of the SM U(l) gauge boson), the final 
decay in the cascade would predominantly be Xi ~~ * lG, 



with two cascades per event, leading to final states with 
77 + E™ lss , where £"™ lss results from the undetected grav- 
itinos. 

Two different classes of gauge-mediated models, de- 
scribed in more detail below, are considered as benchmarks 
to evaluate the reach of this analysis: the minimal GMSB 
model (SPS8) as an example of a complete SUSY model 
with a full particle spectrum and two different variants of 
the GGM model as examples of phenomenological models 
with reduced particle content. 

In the SPS8 model, the only free parameter is the SUSY- 
breaking mass scale A that establishes the nature of the 
observable phenomena exhibited by the low-energy sec- 
tor. The other model parameters are fixed to the follow- 
ing values: the messenger mass M mcss = 2A, the number 
of SU(5) messengers AT 5 = 1, the ratio of the vacuum ex- 
pectation values of the two Higgs doublets tan/? = 15, 
and the Higgs sector mixing parameter fi > 0. The bino 
NLSP is assumed to decay promptly (ctnlsp < 0.1mm). 
For A ~ 200 TeV, the direct production of gaugino pairs 
such as X2 xf or Xi Xi pairs is expected to dominate at 
a LHC centre-of-mass energy of t/s = 7 TeV. The contri- 
bution from gluino and/or squark pairs is below 10% of 
the production cross section due to their high masses. The 
sparticle pair produced in the collision decays via cascades 
into two photons and two gravitinos. Further SM particles 
such as gluons, quarks, leptons and gauge bosons may be 
produced in the cascade decays. The current best limit on 
A in this model is 145 TeV Q. 

Two different configurations of the GGM SUSY model 
are considered in this study, for which the neutralino 
NLSP, chosen to be the bino, and either the gluino or 
the squark masses are treated as free parameters. For the 
squark-bino GGM model all squark masses are treated 
as degenerate except the right-handed up-type squarks 
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whose mass is decoupled (set to inaccessibly large values) . 
For the gluino-bino model all squark masses are decou- 
pled. For both configurations all other sparticle masses are 
also decoupled, leading to a dominant production mode at 
= 7TcV of a pair of squarks in one case and a pair 
of gluinos in the other case. These would decay via short 
cascades into the bino-like neutralino NLSP. Jets may be 
produced in the cascades from the gluino and squark de- 
cays. Further model parameters are hxed to tan /3 — 2 and 
ctnlsp < 0.1 mm. The decay into the wino-like neutralino 
NLSP is possible and was studied by the CMS Collabora- 
tion [Ht 

3. Extra dimensions 

UED models postulate the existence of additional spa- 
tial dimensions in which all SM particles can propagate, 
leading to the existence of a series of excitations for each 
SM particle, known as a Kaluza-Klein (KK) tower. This 
analysis considers the case of a single UED, with compact- 
ihcation radius (size of the extra dimension) R ss 1 TeV -1 . 
At the LHC, the main UED process would be the produc- 
tion via the strong interaction of a pair of first-level KK 
quarks and/or gluons [32j. These would decay via cascades 
involving other KK particles until reaching the lightest KK 
particle (LKP), i.e. the first level KK photon 7*. SM par- 
ticles such as quarks, gluons, leptons and gauge bosons 
may be produced in the cascades. If the UED model is 
embedded in a larger space with N additional eV _1 -sizcd 
dimensions accessible only to gravity [33J, with a (4 + N)- 
dimensional Planck scale (Mb) of a few TeV, the LKP 
would decay gravitationally via 7* — > 7 + G. G represents 
a tower of eV-spaced graviton states, leading to a gravi- 
ton mass between and l/R. With two decay chains per 
event, the final state would contain jj+E™ 
results from the escaping gravitons. Up to l/R 
the branching ratio to the diphoton and E 
close to 100%. As l/R increases, the gravitational decay 
widths become more important for all KK particles and 
the branching ratio into photons decreases, e.g. to 50% 
for 1/R= 1.5 TeV Q. 

The UED model considered here is defined by specifying 
R and A, the ultraviolet cut-off used in the calculation of 
radiative corrections to the KK masses. This analysis sets 
A such that AR = 20 [34j. The 7* mass is insensitive to A, 
while other KK masses typically change by a few per cent 
when varying AR in the range 10 — 30. For 1/R = 1.4 TeV, 
the masses of the first-level KK photon, quark and gluon 
are 1.40 TeV, 1.62 TeV and 1.71 TeV, respectively 



u,lfS where E^ lss 
ITeV, 
™ 1SS final state is 



MRST2007 LO* [40( parton distribution functions (PDFs). 
Signal cross sections were calculated to next-to-leading 
order (NLO) in the strong coupling constant, including 
the resummation of soft gluon emission at next-to-leading- 
logarithmic accuracy (NLO+NLL) • Tne nominal 
cross sections and the uncertainties were taken from an en- 
velope of cross-section predictions using different PDF sets 
and factorisation and renormalisation scales, as described 
in Ref. 46j. In the case of the UED model, cross sections 
were estimated and MC signal samples generated using 
the UED model as implemented at leading order (LO) in 



PYTHIA 6.423 47 



35J with MRST2007 LO* PDFs. 

The "irreducible" background from W(—> Iv) + 77 and 
Z(— > vv) + 77 production was simulated at LO using 
MadGraph 4 Q with the CTEQ6L1 PDFs. Parton 
showering and fragmentation were simulated with PYTHIA. 
NLO cross sections and scale uncertainties were imple- 
mented via multiplicative constants (iVfactors) that re- 
late the NLO and LO cross sections. These have been cal- 
culated for several restricted regions of the overall phase 
space of the Z{— > vv) + 77 and W(— > Iv) + 77 pro- 
cesses 0, HH, and are estimated to be 2.0 ± 0.3 and 3 ± 3 
for the Z{— > vv) + 77 and W{— > Iv) + 77 contributions 
to the signal regions of this analysis, respectively. As de- 
scribed below, all other background sources are estimated 
through the use of control samples derived from data. 

All samples were processed throu gh the GEANT4-based 
simulation of the ATLAS detector [52l IBljj ]. The varia- 
tion of the number of pp interactions per bunch crossing 
(pile-up) as a function of the instantaneous luminosity is 
taken into account by overlaying simulated minimum bias 
events according to the observed distribution of the num- 
ber of pile-up interactions in data, with an average of ~ 10 
interactions. 

5. ATLAS detector 

The ATLAS detector [54| is a multi-purpose apparatus 
with a forward-backward symmetric cylindrical geometry 
and nearly Att solid angle coverage. Closest to the beam- 
line are tracking devices comprising layers of silicon-based 
pixel and strip detectors covering I77I < 2.f0] and straw- 
tube detectors covering |?7| < 2.0, located inside a thin 
superconducting solenoid that provides a 2T magnetic 
field. Outside the solenoid, fine-granularity lead/liquid- 
argon electromagnetic (EM) calorimeters provide cover- 
age for |?7| < 3.2 to measure the energy and position of 
electrons and photons. A presampler, covering \r/\ < 1.8, 
is used to correct for energy lost upstream of the EM 



4. Simulated samples 

For the GGM model, the SUSY mass spectra were cal- 
culated using SUSPECT 2.41 [36| and SDECAY 1.3 37]; for 
the SPS8 model, the SUSY mass spectra were calculated 
using ISA JET 7.80 The Monte Carlo (MC) SUSY sig- 
nal samples were produced using Herwig++ 2.5.1 39] with 



1 ATLAS uses a right-handed coordinate system with its origin 
at the nominal interaction point (IP) in the centre of the detector 
and the z-axis along the beam pipe. The x-axis points from the 
IP to the centre of the LHC ring, and the j/-axis points upward. 
Cylindrical coordinates (r, <j>) are used in the transverse plane, (j> 
being the azimuthal angle around the beam pipe. The pseudorapidity 
is defined in terms of the polar angle 9 as rj = — lntan(6>/2). 
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calorimeter. An iron/scintillating-tile hadronic calorime- 
ter covers the region \r)\ < 1.7, while a copper /liquid- argon 
medium is used for hadronic calorimeters in the end-cap re- 
gion 1.5 < \?]\ < 3.2. In the forward region 3.2 < \ri\ < 4.9 
liquid-argon calorimeters with copper and tungsten ab- 
sorbers measure the electromagnetic and hadronic energy. 
A muon spectrometer consisting of three superconduct- 
ing toroidal magnet systems each comprising eight toroidal 
coils, tracking chambers, and detectors for triggering sur- 
rounds the calorimeter system. 

6. Reconstruction of candidates and observables 

The reconstruction of converted and unconverted pho- 
tons and of electrons is described in Refs. [55[ and (Hi}, re- 
spectively. Photon candidates were required to be within 
|?7| < 1.81, and to be outside the transition region 1.37 < 
|?7| < 1.52 between the barrel and end-cap calorimeters. 
Identified on the basis of the characteristics of the lon- 
gitudinal and transverse shower development in the EM 
calorimeter, the analysis made use of both "loose" and 
"tight" photons [EBj]. In the case that an EM calorimeter 
deposition was identified as both a photon and an electron, 
the photon candidate was discarded and the electron can- 
didate retained. In addition, converted photons were re- 
classified as electrons if one or more candidate conversion 
tracks included at least one hit from the pixel layers. Giv- 
ing preference to the electron selection in this way reduced 
the electron-to-photon fake rate by 50-60% (depending on 
the value of rj) relative to that of the prior lfb -1 analy- 
sis [l[, while preserving over 70% of the signal efficiency. 
Finally, an "isolation" requirement was imposed. After 
correcting for contributions from pile-up and the deposi- 
tion ascribed to the photon itself, photon candidates were 
removed if more than 5 GeV of transverse energy was ob- 
served in a cone of J (A77) 2 + (A</>) 2 < 0.2 surrounding 
the energy deposition in the calorimeter associated with 
the photon. 

The measurement of the two-dimensional transverse 
momentum vector p™ ss (and its magnitude E™ ss ) was 
based on energy deposits in calorimeter cells inside three- 
dimensional clusters with |?7| < 4.9 and was corrected for 
contributions from muons, if any [57] . The cluster en- 
ergy was calibrated to correct for the different response 
to electromagnetically- and hadronically-induced showers, 
energy loss in dead material, and out-of-cluster energy. 
The contribution from identified muons was accounted for 
by adding in the energy derived from the properties of 
reconstructed muon tracks. 

Jets were reconstructed using the anti-fct jet algo- 
rithm [58| with radius parameter R = 0.4. They were 
required to have pt > 20 GeV and | ^7 1 < 2.8 [59] . 

Two additional observables of use in discriminating SM 
backgrounds from potential GMSB and UED signals were 
defined. The total visible transverse energy Ht was cal- 
culated as the sum of the magnitude of the transverse mo- 
menta of the two selected photons and any additional lep- 



tons and jets in the event. The photon- E T 11SS separation 
A0(7, E T USS ) was defined as the azimuthal angle between 
the missing transverse momentum vector and either of the 
two selected photons, with A0 m j n (7, _E™ 1SS ) the minimum 
value of A</>(7, E™ lss ) of the two selected photons. 

7. Data analysis 

The data sample, corresponding to an integrated lumi- 
nosity of (4.8 ± 0.2) fb -1 [6fj, |6l|, was selected by a trigger 
requiring two loose photon candidates with Et > 20 GeV. 
To ensure the event resulted from a beam collision, events 
were required to have at least one vertex with five or more 
associated tracks. Events were then required to contain 
at least two tight photon candidates with Et > 50 GeV, 
which MC studies suggested would provide the greatest 
separation between signal and SM background for a broad 
range of the parameter space of the new physics scenar- 
ios under consideration in this search. A total of 10455 
isolated 77 candidate events passing these selection re- 
quirements were observed in the data sample. The Et 
distribution^] of the leading and sub-leading photon for 
events in this sample are shown in Figs. [T] and [2j Also 
shown are the Et spectra obtained from GGM MC sam- 
ples for m~ g = 1000 GeV and m^o = 450 GeV, from SPS8 
MC samples with A = 190 TeV, and from UED MC sam- 
ples for 1/R = 1.3 TeV, representing model parameters 
near the expected exclusion limit. Figures[3]and[4]show the 
Ht and A0 m i n (7, E™ lss ) distributions of selected diphoton 
events, with those of the same signal models overlaid. 

To maximise the sensitivity of this analysis over a wide 
range of model parameters that may lead to different kine- 
matic properties, three different signal regions (SRs) were 
defined based on the observed values of E!f lss , Ht and 
A0 m i n (7, i?™ lss ). SR A, optimised for gluino/squark pro- 
duction with a subsequent decay to a high-mass bino, re- 
quires large E™ 1SS and moderate Ht- SR B, optimised 
for gluino/squark production with a subsequent decay 
to a low-mass bino, requires moderate E™ lss and large 
Ht- SR C, optimised for the electroweak production of 
intermediate-mass gaugino pairs that dominates the SPS8 
cross section in this regime, requires moderate E™ lss but 
makes no requirement on Ht- In addition, a requirement 
of A(/) m j n (7, E™ ms ) > 0.5 was imposed on events in SR A 
and C; for the low-mass bino targeted by SR B, the sepa- 
ration between the photon and gravitino daughters of the 
bino is too slight to allow for the efficient separation of 



2 An excess of events relative to a smoothly-falling distribution 
of the leading-photon spectrum was observed for Et ~ 285 GeV. 
Searching over the range 100 GeV < Et < 500 GeV, a significance 
of 1.9 c was found using BumpHunter lO'J . while the local signifi- 
cance was found to be 3.1 a. No correlation between the excess and 
the LHC running period or luminosity was observed. A comparison 
of other observables (e.g. diphoton mass, _E™ 1SS , leading-photon r/, 
A</>(7i,72)) between the excess and sideband regions exhibited no 
appreciable differences. It was concluded that the observed excess of 
events is compatible with a statistical fluctuation. 
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signal from background through the use of this observ- 
able. The selection requirements of the three SRs are sum- 
marised in Table [TJ Of the three SRs, SR A provides the 
greatest sensitivity to the UED model, and is thus the SR 
used to test this model. 



Table 1: Definition of the three SRs (A, B and C) based on the 
quantities E™ iss , H T and A?J min (7, E™ iss ). 





SR A 


SR B 


SR C 


E™ lss > 


200 GeV 


100 GeV 


125 GeV 


H T > 


600 GeV 


1100 GeV 




A0 min ( 7 ,^ iss ) > 


0.5 




0.5 



Table [5] shows the numbers of events remaining after 
several stages of the selection. A total of 117, 9 and 7293 
candidate events were observed to pass all but the E^ lss 
requirement of SR A, B and C, respectively. After im- 
posing the final E™ lss requirement, no events remained for 
SR A and B, while two events remained for SR C. 



Table 2: Samples of selected events at progressive stages of the se- 
lection. Where no number is shown the cut was not applied. 



Triggered events 


1166060 


Diphoton selection 


10455 




A 


B 


C 


A0 min (7,ST lss ) requirement 


7293 




7293 


Ht requirement 


117 


9 




E™ lss requirement 
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Figure [5] shows the E™ lss distribution for SR C, the 
expected contributions from the SPS8 MC sample with 
A = 190 TeV, and estimated background contributions 
from various sources (described below). 

8. Background estimation 

Following the procedure described in Ref. [63[, the con- 
tribution to the large E™ lss diphoton sample from SM 
sources can be grouped into three primary components. 
The first of these, referred to as "QCD background" , arises 
from a mixture of processes that include 77 production as 
well as 7 + jet and multijet events with at least one jet 
mis-reconstructed as a photon. The second background 
component, referred to as "EW background", is due to 
W + X and tt events (here "A" can be any number of 
photons or jets), and where mis- reconstructed photons can 
arise from electrons and jets, and for which final-state neu- 
trinos produce significant E™ lss . The QCD and EW back- 
grounds were estimated via dedicated control samples of 
data events. The third background component, referred 
to as "irreducible" , consists of W and Z bosons produced 
in association with two real photons, with a subsequent 
decay into one or more neutrinos. 

To estimate the QCD background from 77, 7 + jet, and 
multijet events, a "QCD control sample" was selected from 
the diphoton trigger sample by selecting events for which 




Data 2011 (^i = 7TeV) 
GGM m(i) = 1000 GeV x1 

m(x°) = 450 GeV 
SPS8A=190TeVx100 

UED1/R=1.3TeVx100 
Ldt = 4.8 fb 1 
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Figure 1: The &p spectrum of the leading photon in the 77 candi- 
date events in the data (points, statistical uncertainty only) together 
with the spectra from simulated GGM (nig = 1000 GeV, m-o = 
450 GeV), SPS8 (A = 190 TeV), and UED [X/R = 1.3 TeV) samples 
after the diphoton requirement. The signal samples are scaled by a 
factor of 100 for clarity. 



at least one of the photon candidates passes the loose but 
not the tight photon identification. Events with electrons 
were vetoed to remove contamination from W — > ev de- 
cays. The Ht and A</>( 7 , E™ lss ) requirements associated 
with each of the three SRs were then applied, yielding 
three separate QCD samples, or "templates" . An estimate 
of the QCD background contamination in each SR was 
obtained from imposing the £"™ lss requirement associated 
with the given SR upon the corresponding QCD template, 
after normalising each template to the diphoton data with 
E™ isB < 20 GeV from the given SR. This yielded a QCD 
background expectation of 0.85 ± 0.30(stat) events for SR 
C. No events above the corresponding E™ lss requirement 
were observed for the A and B control samples, yielding 
an estimate of events with a 90% confidence-level (CL) 
upper limit of less than 1.01 and 1.15 background events 
for SR A and SR B, respectively. 

To improve the constraint on the estimated background 
for SRs A and B, a complementary method making use of 
Ht sidebands of the QCD control sample was employed. 
The Ht requirement applied to the QCD templates of SR 
A and B was relaxed in three steps: to 400 GeV, 200 GeV 
and GeV for the SR A control sample, and to 800 GeV, 
400 GeV and 200 GeV for the SR B control sample. For 
each SR, the E!f lss distribution of each of these relaxed 
control samples was scaled to the diphoton E™ lss distri- 
bution for S™ iss < 20 GeV of the given SR, yielding a 
series of three expected values for the QCD background 
as a function of the applied Ht requirement. The com- 
plementary estimate for the background contribution to 
the signal region employed a parabolic extrapolation to 
the actual Ht requirement used for the analysis (600 GeV 
and 1100 GeV for SRs A and B, respectively); a linear 
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GGM m(i) = 1000GeV x1 

m(jp = 450 GeV 
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UED1/R=1.3TeVx100 
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Figure 2: The Ex spectrum of the sub-leading photon in the 
77 candidate events in the data (points, statistical uncertainty 
only) together with the spectra from simulated GGM (mg = 
1000 GeV, m-o = 450 GeV), SPSS (A = 190TeV), and UED (l/R = 
1.3 TeV) samples after the diphoton requirement. The signal samples 
are scaled by a factor of 100 for clarity. 



Figure 3: The Hx spectrum of 77 candidate events in the data 
(points, statistical uncertainty only) together with the spectra from 
simulated GGM (m s = 1000 GeV, m- = 450 GeV), SPS8 (A = 
190TeV), and UED (l/R = 1.3TeV) samples after the diphoton 
requirement. The signal samples are scaled by a factor of 100 for 
clarity. 



fit yielded a significantly lower background estimate for 
both SRs. The parabolic extrapolation yielded conserva- 
tive upper estimates of 0.14 and 0.54 events for SRs A and 
B, respectively. The overall QCD background estimates 
for SRs A and B were taken to be 0.07 ± 0.07(syst) and 
0.27 ± 0.27(syst) events, respectively, half of the value of 
this upper estimate, with systematic uncertainty assigned 
to cover the entire range between and the upper esti- 
mate. The choice of a parabolic function constrained by 
three Hx points does not permit an estimation of statisti- 
cal uncertainty on the extrapolation. 

Other sources of systematic uncertainty in the estimated 
QCD background were considered. Using the E™ lss dis- 
tribution from a sample of Z — > e + e~ events instead of 
that of the QCD sample yielded estimates of 0, and 
0.15 events for the SRs A, B and C, respectively. The 
difference between this estimate and that of the QCD 
sample was incorporated as a systematic uncertainty of 
±0.71 on the SR C QCD background estimate. Making 
use of the alternative ranges 5 GeV < i?™ lss < 25 GeV 
and 10 GeV < E™ iss < 30 GeV over which the QCD sam- 
ple was normalized to the 77 sample resulted in a further 
systematic uncertainty of ±0.03 events on the QCD back- 
ground estimate for SR C. The resulting QCD background 
estimates for the three SRs, along with their uncertainties, 
are compiled in Table [3] 

The EW background, from W + X and ti events, was 
estimated via an "electron-photon" control sample com- 
posed of events with at least one tight photon and one 
electron, each with Ex > 50 GeV, and scaled by the 
probability for an electron to be mis-reconstructed as a 
tight photon, as estimated from a "tag-and-probe" study 
of the Z boson in the ee and ej sample. The scaling 



factor varies between 2.5% (0 < \r)\ < 0.6) and 7.0% 
(1.52 < < 1.81), since it depends on the amount of 
material in front of the calorimeter. Events with two or 
more tight photons were vetoed from the control sample 
to preserve its orthogonality to the signal sample. In case 
of more than one electron, the one with the highest px was 
used. 

After applying corresponding selection requirements on 
H T , A</>(7, ££ liss ) and £^ iss , a total of 1, 3 and 26 
electron-photon events were observed for SRs A, B and 
C, respectively. After multiplying by the ^-dependent 
electron-to-photon mis-reconstruction probability, the re- 
sulting EW background contamination was estimated to 
be 0.03 ± 0.03, 0.09 ± 0.05 and 0.80 ± 0.16 events for SRs 
A, B and C, respectively, where the uncertainties are sta- 
tistical only. 

The systematic uncertainty on the determination of the 
electron-to-photon mis-reconstruction probability is as- 
sessed by performing an independent tag-and-probe anal- 
ysis with looser electron Ex and identification require- 
ments. Differences with the nominal tag-and-probe anal- 
ysis are taken as systematic uncertainty on the EW back- 
ground estimate, resulting in relative systematic uncer- 
tainties of ±6.9%, ±7.1% and ±10.0% for SRs A, B and 
C, respectively. MC studies suggest that approximately 
25% of the EW background involves no electron-to-photon 
mis-reconstruction, and thus are not accounted for with 
the electron-photon control sample. These events, how- 
ever, typically involve a jet-to-photon mis-reconstruction 
(for example, an event with one radiated photon and 
a hadronic r decay with an energetic leading tt mis- 
reconstructed as a photon), and are thus potentially ac- 
counted for in the QCD background estimate. A relative 
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Table 3: The expected number of 77 events for each of the three signal regions. The uncertainties are statistical, arising from the limited 
numbers of events in the control samples, and systematic, the details of which are given in the text. For the irreducible background, the 
statistical uncertainty is due to limited numbers of events in the corresponding MC samples. 





SR A 


SR B 


SR C 


QCD 

Electroweak 

iv) + 77 
Z{— > vv) + 77 


0.07 ±0.00 ±0.07 
0.03 ±0.03 ±0.01 

< 0.01 

< 0.01 


0.27 ±0.00 ±0.27 
0.09 ±0.05 ±0.02 

< 0.01 

< 0.01 


0.85 ±0.30 ±0.71 
0.80 ±0.16 ±0.22 
0.18 ±0.13 ±0.18 
0.27 ±0.09 ±0.04 


Total 


0.10 ±0.03 ±0.07 


0.36 ±0.05 ±0.27 


2.11 ±0.37 ±0.77 


Observed events 








2 




Figure 4: The minimum A<^>(7, £?™ lss ) spectrum of 77 candidate 
events in the data (points, statistical uncertainty only) together 
with the spectra from simulated GGM (nig = 1000 GeV, m^o = 

450 GeV), SPS8 (A = 190 TeV), and UED (1/R = 1.3 TeV) samples 
after the diphoton requirement. The signal samples are scaled by a 
factor of 100 for clarity. 



Figure 5: Ej? lss spectra in SR C for the 77 candidate events in data 
(points, statistical uncertainty only) and the estimated QCD back- 
ground (normalised to the number of 77 candidates with _E™ 1SS < 
20 GeV), the W(— > ev) + jets/7 an d it(— > ev) + jets backgrounds 
as estimated from the electron-photon control sample, and the ir- 
reducible background of Z(— > vv) + 77 and W(— > iv) + 77. The 
hatched region represents the extent of the uncertainty on the total 
background prediction. Also shown is the expected signal from the 
SPS8 (A = 190 TeV) sample. 



systematic uncertainty of ±25% is conservatively assigned 
to the EW background estimates for all three SRs to ac- 
count for this ambiguity. The resulting EW background 
estimates for the three SRs, along with their uncertainties, 
are compiled in Table [21 

The contribution of the irreducible background from the 
Z(—> w)+77 and W(— > £v)+jj processes was estimated 
using MC samples. It was found to be negligible for SRs 
A and B, and estimated to be 0.46 ± 0.16 ± 0.19 events 
for SR C, where the first uncertainty is due to the limited 
number of events in the MC sample and the second to 
the uncertainty on the applied iVfactor. These estimates, 
along with the resulting estimates for the total background 
from all sources, are reported in Table [31 

The contamination from cosmic-ray muons, estimated 
using events triggered in empty LHC bunches, was found 
to be negligible. 



9. Signal efficiencies and systematic uncertainties 

Signal efficiencies were estimated using MC simulation. 
GGM signal efficiencies were estimated over an area of 
the GGM parameter space that ranges from 800 GeV to 
1300 GeV for the gluino or squark mass, and from 50 GeV 
to within 10 GeV of the gluino or squark mass for the neu- 
tralino mass. For SR A the efficiency increases smoothly 
from 1.2% to 25% for (m^m^o) = (800, 50) GeV to 
(1300, 1280) GeV, but then drops to 20% for the case for 
which the gluino and neutralino masses are only separated 
by 10 GeV. For SR B the efficiency increases smoothly 
from 2.8% to 26% for (m s ,m^o) = (800, 790) GeV to 
(1300, 50) GeV. The SPS8 signal efficiency in SR C in- 
creases smoothly from 5.9% (A = 100 TeV) to 21% 
(A = 250 TeV). For SR A the UED signal efficiency 
increases smoothly from 28% (1/R = 1.0 TeV) to 37% 
(1/R = 1.5 TeV). 

The various relative systematic uncertainties on the 
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GGM, SPS8 and UED signal cross sections are summarised 
in Table [4] for the chosen reference points: (mj,m^) = 
(1000, 450) GeV for GGM, A = 190 TeV for SPSs/and 
1/R = 1.3 TeV for UED. The uncertainty on the luminos- 
ity is ±3.9 % [go|, 61 1 . The efficiency of the required dipho- 
ton trigger was estimated using a single photon trigger, 
the efficiency of which was determined using a bootstrap 
method The result is 99. 8^:8% for events passing 

the diphoton selection. To estimate the systematic un- 
certainty due to the unknown composition of the data 
sample, the trigger efficiency was also evaluated on MC 
events using mis-reconstructed photons from filtered mul- 
tijet samples and photons from signal (GGM, SPS8 and 
UED) samples. A conservative systematic uncertainty of 
±0.5 % was derived from the difference between the ob- 
tained efficiencies. Uncertainties on the photon selection, 
the photon energy scale, and the detailed material com- 
position of the detector, as described in Ref. [63j], result 
in an uncertainty of ±4.4 % for the GGM, SPS8 and UED 
signals. The uncertainty due to the photon isolation re- 
quirement was estimated by varying the energy leakage 
and the pile-up corrections independently, resulting in an 
uncertainty of ±0.9%, ±0.2% and ±0.4% for the GGM, 
SPS8 and UED signals, respectively. The influence of pile- 
up on the signal efficiency, evaluated by scaling the number 
of pile-up events in the MC simulation by a factor of 0.9 
(chosen to reflect the range of uncertainty inherent in es- 
timating and modeling the effects of pile-up) leads to a 
systematic uncertainty of ±0.8% (GGM), ±0.5% (SPS8) 
and ±0.5% (UED). Systematic uncertainties due to the 
jrjmiss reconstruction, estimated by varying the cluster en- 
ergies and the i?™ ISS resolution between the measured per- 
formance and MC expectations [57) . contribute an un- 
certainty of ±0.1/0.5% to ±5.3/16.1% (GGM, SR A/B), 



Table 4: Relative systematic uncertainties on the expected signal 
yield for the GGM model with m- g = 1000 GeVandm -o = 450 GeV, 
the SPS8 model with A = 190 TeV, and the UED model with l/R = 
1.3 TeV. For the GGM model, when the uncertainty differs for SRs A 
and B, it is presented as SRA/SRB. No PDF and scale uncertainties 
are given for the UED case as the cross section is evaluated only to 
LO. 



Source of uncertainty 


Uncertainty 






GGM 


SPS8 


UED 


Integrated luminosity 


3.9% 


3.9% 


3.9% 


Trigger 


0.5% 


0.5% 


0.5% 


Photon identification 


4.4% 


4.4% 


4.4% 


Photon isolation 


0.9% 


0.2% 


0.4% 


Pile-up 


0.8% 


0.5% 


0.5% 


^miss reconduction 


3.9/1.1% 


2.8% 


1.5% 




0.0/2.1% 




0.4% 


Signal MC sample size 


3.0% 


2.1% 


1.4% 


Total signal uncertainty 


7.6/7.1% 


6.8% 


6.3% 


PDF and scale 


31% 


5.5% 




Total 


32% 


8.7% 


6.3% 
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Figure 6: Expected and observed 95% CL lower limits on the gluino 
mass as a function of the neutralino mass in the GGM model with 
a bino-like lightest neutralino NLSP (the grey area indicates the 
region for which the gluino mass is less than the bino mass, which is 
not considered here). The other sparticle masses are assumed to be 
decoupled. Further model parameters are tan/3 = 2 and ctnlsp < 
0.1mm. The previous ATLAS limit Q is also shown. 



±1.6% to ±9.7% (SPS8) and ±0.9% to ±2% (UED). Sys- 
tematic uncertainties due to the ifr reconstruction, esti- 
mated by varying the energy scale and resolution of the 
individual objects entering Ht, are below ±0.3% (GGM, 
SR A), ±0.1% to ±7.3% (GGM, SR B) and ±1.1% to 
±0.1% (UED). The systematic uncertainties from E!f lss 
and Ht are taken to be fully correlated. Added in 
quadrature, the total systematic uncertainty on the sig- 
nal yield varies between ±6% and ±20% (GGM), ±6% 
and ±15% (SPS8), and ±6% and ±7% (UED). 

The PDF and factorisation and renormalisation scale 
uncertainties on the GGM (SPS8) cross sections were eval- 
uated as described in Section 4, leading to a combined 
systematic uncertainty between ±23-39%, ±29-49% and 
±4.7-6.4% for the GGM (gluino), GGM (squark) and 
SPS8 models, respectively. The different impact of the 
PDF and scale uncertainties on the GGM and SPS8 yields 
is related to the different production mechanisms in the 
two models (see Section [2]). In the case of UED, the 
PDF uncertainties were evaluated by using the MSTW2008 
L0 65j PDF error sets in the LO cross-section calculation 
and are about ±4%. The scale of a s in the LO cross sec- 
tion calculation was increased and decreased by a factor 
of two, leading to a systematic uncertainty of ±4.5% and 
±9 %, respectively. NLO calculations are not yet available, 
so the LO cross sections were used for the limit calcula- 
tion without any theoretical uncertainty, and the effect 
of PDF and scale uncertainties on the final limit is given 
separately. 
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Figure 7: Expected and observed 95 % CL lower limits on the squark 
mass as a function of the neutralino mass in the GGM model with 
a bino-like lightest neutralino NLSP (the grey area indicates the 
region for which the squark mass is less than the bino mass, which is 
not considered here). The other sparticle masses are assumed to be 
decoupled. Further model parameters are tan/3 = 2 and ctjmlsp < 
0.1 mm. 



10. Results 

No evidence for physics beyond the SM was observed in 
any of the SRs. Based on the numbers of observed events 
in SR A, B and C and the background expectation shown 
in Table [3J 95 % CL upper limits are set on the numbers 
of events in the different SRs from any scenario of physics 
beyond the SM using the profile likelihood and CL S pre- 
scriptions [66[| . Uncertainties on the background and signal 
expectations are treated as Gaussian-distributed nuisance 
parameters in the maximum likelihood fit, resulting in ob- 
served upper limits of 3.1, 3.1 and 4.9 events for SRs A, 
B and C, respectively. These limits translate into 95% 
upper limits on the visible cross section for new physics, 
defined by the product of cross section, branching ratio, 
acceptance and efficiency for the different SR definitions, 
of 0.6, 0.6 and 1.0 fb, respectively. Because the observed 
numbers of events are close to the expected numbers of 
background events for all three SRs, expected limits on 
the numbers of events from and visible cross section for 
new physics are, to the quoted accuracy, identical to the 
observed limits. 

Limits are also set on the GGM squark and gluino 
masses as a function of the bino-like neutralino mass, mak- 
ing use of the SR (A or B) that provides the most stringent 
expected limit for the given neutralino mass. Figures [6] 
and [7] show the expected and observed lower limits on the 
GGM gluino and squark masses, respectively, as a function 
of the neutralino mass. Three observed-limit contours are 
shown, corresponding to the nominal assumption for the 
SUSY production cross section as well as those derived 
by reducing and increasing the cross section by one stan- 
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Figure 8: Expected and observed 95 % CL upper limits on the spar- 
ticle production cross section in the SPS8 model, and the NLO cross- 
section prediction, as a function of A and the lightest neutralino and 
chargino masses. Further SPS8 model parameters are M meS s = 2A, 
N$ = 1, tan/3 = 15 and ctnlsp < 0.1mm. Limits are set based on 
SR C. 



dard deviation of theoretical uncertainty (the combined 
uncertainty due to the PDFs and renormalisation and fac- 
torisation scales). For comparison the lower limits on the 
GGM gluino mass from ATLAS based on 1 ftr 1 from 
2011 are also shown. 

Including all sources of uncertainty other than the the- 
oretical uncertainty, 95 % CL upper limits on the cross 
section of the SPS8 model are derived from the SR C 
result and displayed in Fig. [8] for the range A = 100- 
250 TeV along with the overall production cross section 
and its theoretical uncertainty. For illustration the cross- 
section dependence as a function of the lightest neutralino 
and chargino masses is also shown. 

Figure [5] shows the limit on the cross section times 
branching ratio for the UED model as a function of the 
compactification scale 1 /R, derived from the result of SR 
A. A 95% CL lower limit of 1/R > 1.40 TeV is set. For 
illustration the cross-section dependence as a function of 
the KK quark and KK gluon masses is also shown. Again, 
neither PDF nor scale uncertainties are included when cal- 
culating the limits; including PDF and scale uncertainties, 
computed at LO, in the limit calculation degrades the limit 
on 1/R by a few GeV. 

11. Conclusions 



A search for events with two photons and substantial 
performed using 4.8 fb -1 of 7 TeV pp collision data 
recorded with the ATLAS detector at the LHC, is pre- 
sented. The sensitivity to different new physics mod- 
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Figure 9: Expected and observed 95 % CL upper limits on the KK 
particle production cross section times branching ratio to two pho- 
tons in the UED model, and the LO cross-section prediction times 
branching ratio, as a function of l/R and the KK quark (Q*) and 
KK gluon (g*) masses. The ±1<t expected-limit error band overlaps 
the observed limit contour and is too narrow to be distinguished. 
No error is shown for the UED cross section since the cross-section 
calculation is available only to LO (see text for further discussion). 
The UED model parameters are TV = 6, M D = 5 TeV and AR = 20. 
Limits arc set based on SR A. 



els producing this final state was optimised by defining 
three different SRs. No significant excess above the ex- 
pected background is found in any SR. The results are 
used to set model-independent 95 % CL upper limits on 
possible contributions from new physics. In addition, 
under the GGM hypothesis, considering cross sections 
one standard deviation of theoretical uncertainty below 
the nominal value, a lower limit on the gluino/squark 
mass of 1.07TeV/0.87TeV is determined for bino masses 
above 50 GeV. Under similar assumptions, a lower limit of 
196 TeV is set on the SUSY-breaking scale scale A of the 
SPS8 model. Considering nominal values of the leading- 
order UED cross section, a lower limit of 1.40 TeV is set on 
the UED compactification scale l/R. These results pro- 
vide the most stringent tests of these models to date. 
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